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Abstract. The structure and dynamics of an n-particle system are described with 
coupled nonlinear Heisenberg's commutator equations where the nonlinear terms 
are generated by the two-body interaction that excites the reference vacuum via 
particle-particle and particle-hole excitations. Nonperturbative solutions of the 
system are obtained with the use of dynamic linearization approximation and cluster 
transformation coefficients. The dynamic linearization approximation converts the 
commutator chain into an eigenvalue problem. The cluster coefficients factorize the 
matrix elements of the ('nj-particles or particle- hole systems in terms of the matrix 
elements of the (n- 1 )-systems coupled to a particle-particle, particle-hole, and hole- 
hole boson. Group properties of the particle-particle, particle-hole, and hole-hole 
permutation groups simplify the calculation of these coefficients. The particle-particle 
vacuum-excitations generate superconductive diagrams in the dynamics of 5-quarks 
systems. Applications of the model to fermionic and bosonic systems are discussed. 
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1. Introduction 

In the Heisenberg's picture the time evolution of a system of particles is described 
by a commutator equation which involves the n-body Hamiltonian and the creation 
operators of the n-body ground- and excited-modes. However, the excitations of 
the reference vacuum, resulting from the scattering of particles from the vacuum 
to higher states (the particle-hole and particle-particle excitations), are completely 
neglected in this formulation, namely, the time evolution of the n-body-modes 
is described by a linearized Equation of Motion (EOM) [1] which involves only 
valence particles. Recently the original Heisenberg formulation of the n-body- 
dynamics has been generalized within the framework of the Dynamic-Correlation 
Model (DCM) [2] for fermions, the Boson Dynamic-Correlation Model (BDCM) [3] for 
bosons, and the Superconductive Dynamic- Correlation Model (SDCM) which describes 
superconductive- and polarization-effects [4]. In these models, the inclusion of the 
structure and the dynamics of odd/even quasi-particles into the calculations of the 
excitations of the model-vacuum led to the modification of the formal strcture of the 
original Heisenberg's commutator equation. The new dynamics system is characterized 
by a set of coupled commutator equations which involve simultaneously the excitations of 
the valence quasi-particles and the excitations of the Intrinsic-Supcrcunductive- Vacuum 
States (ISVSs) i.e.: mixed-mode states formed by coupling the valence quasi-particles 
to the excitations of the vacuum (particle-particle and particle-hole). As in Ref. [4], 
the resulting mixed-mode states are classified in terms of the following mixed-mode 
wave functions: (a) valence particles coupled to ISVSs formed by include particle- 
hole, particle-particle and hole-hole vacuum-excitation modes that have the same parity 
of the valence particles; (b) valence particles and ISVSs formed by exciting particle- 
hole, particle-particle and hole-hole having a parity different from that of the valence 
particles. The superconductive vacuum states (b) arc not considered by perturbative 
theories, although they may be associated to the creation of virtual particles [2, 3] 
giving important contribution to dynamic theories. The superconductive vacuum 
excitation modes are important mainly at high densities [4], although the overlap of 
the particle-particle excitation modes may infiuence, due to the strong Pauli blocking 
effects in the ISVSs, the particle-hole excitation mechanism also at medium energies. 
In this paper we disregard the superconductive effects, which will be the subject of 
a paper in preparation, and discuss only the effect of the particle-hole excitation 
mechanism on ground state properties of atomic, nuclear, and quark systems. In 
principle, the system of commutator equations may be solved perturbativcly by means 
of substitution method which consists in inserting the higher order commutators into 
the previous calculated commutators. However, such perturbative method is usually 
not convergent for systems of strong interacting particles and, therefore, will not be 
discussed here. We will study nonperturbative solutions based on the following: (a) 
a linearization ansatz motivated by the consideration that in the low energy domain 
only few particle excitations are energetically representative while the higher excitation 
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modes (with admixture probability less than 0.1 %) may be omitted from the model 
space; (b) cluster transformation coefficients which provide an exact solution for the 
complex n-body matrix elements that are input of the commutator equations. The 
wave functions, solutions of the calculated collective eigenvalue equation are classified 
as in Refs. [2, 3] in terms of Configuration Mixing Wave Functions (CMWFs). In 
Section II, we calculate the commutator chain for one particle-, two particle-, and three 
particle- creation operators. The generalization of the model to hole- and (particle-hole)- 
valence operators is not given explicitly because of easy extrapolation from the present 
particle formulation. The commutator chain is then linearized and converted into an 
eigenvalue problem which is solved calculating the n-body matrix elements within the 
cluster transformation coefficients. In Section III, solutions of the nonlinear eigenvalue 
equation for one valence particle/hole are obtained and the HPS constants of hydrogen- 
and lithium-like heavy ions are calculated. In Section IV, the commutator chain is solved 
for two and three valence particles and the resulting CMWFs are used to calculate the 
matter distribution of halo nuclei. In Section V, an application to the three dressed quark 
systems is discussed in order to provide new theoretical data for the polarizability of 
the proton. For a measurement of the proton polarization, novel ultra/intense lasers 
may be useful, like the PHELIX-petawatt laser presently under construction at GSI, 
have to be used. A direct way for a precision measurement seems to be in reach by the 
latest improvements of energetic multi-MeV photon sources using laser backscattering 
at electron storage rings. 



2. The nonlinear commutator chain 



The nonlinear commutator chain is an extension of the Heisenberg's equation. It 
allows us to address the situation where valence clusters and core clusters are almost 
energetically degenerated and may, therefore, coexist. We introduce valence systems 
formed of either neutrons or protons, and define the valence states as 

|r,,({n} particles)) = ^ xj^Jj^ N^^'Jj^ 4(a„J.; J)|0), (1) 

(n) 

where n = 1, 2, 3 and J is the total spin. The N'^; 

rJn normalization constants, 

■^a^Jjn mode amplitudes, and cti denotes the other quantum numbers of the valence 
particles. The creation operator is defined by 

A\{ain; J) = a]^ (2) 
for one valence particle (fermion); 

4 (as J; J) = (4, ® a]y (3) 
for one valence pair (boson); 

4(«3^i; J) = ((4, ® a]y^ ® al)' (4) 

for three valence particles (fermion). Hence, in Eq. (|l]), J„ = ji for one particle, Jn = J 
for two particles, and J„ = Ji for three particles. 
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In the literature, various approaches, as reported in Ref. [3], have been proposed for 
including core excitations. In Shell-Model calculations the residual interaction between 
the valence particles excites the valence pairs to higher single particle states, leaving 
the vacuum in the J = ground state. In Ref. [5] the core excitations have been 
included in the Shell-Model calculations through introducing the coupling between the 
valence particles and the collective J = 2^ core-excited states. In the DCM and BDCM 
the core excitation is included through coupling the valence fermionic/bosonic states, 
Eqs. (H, H) to intrinsic bosonic states corresponding to particle-hole excitations of the 
nuclear core. In this paper, we consider the following mixed-mode fermionic/bosonic 
states: a) valence fermionic/bosonic states coupled to the dynamic particle-hole states 
of normal parity; b) valence fermionic/bosonic states coupled to the dynamic particle- 
hole states of non-normal parity. The particle-hole coupling is implemented through the 
two-body force V , such that 

if = ^ aj^ a„ + ^ ^ v^p^^s a{; ^5 a^ = Hq + V, (5) 

where Va^^^, are the matrix elements of the realistic two-body potential which has two 

parts: a central part and a tensor part. The tensor component of the realistic two- 
body potential shapes the many-body Hamiltonian and in particular its long tail acting 
between the valence fermions/bosons and causing simultaneously the excitation of the 
particles to high shell-model states and the deformation of the nuclear core. 

The mixed-mode {1 — 2 — 3} partides{l' — 2'} bosons states, denoted | Fj), can 
therefore be expanded as follows: 

\T'ji{n = 1-2-3 partzdes} - {l' - 2'} bosons)) = 4"jjX:!j„4(«n Jn; J) 

+ Yl ^i!ti^,ij„+i,^i!ti,J„J„+i,^n+l'(an+l'(^^+l'; J)) 

where 1,2,3 labels the particles and l',2' label the bosons (particle-hole pairs). The 
operators 

4+l,(a„+r(JlJ„+lO; J) = ([4(«n(Jn; J) ® {al^, ® (7) 
A\^2'{(^n+2'{JnJn+vJn+2')] J) = ( [AJ,, («„ ( J„; J))® 

® bl^j'-^" ® («k+3 ® 6;„^J^"+2')''<"+^' "+^'']-'"'^^"+^' "+^'0'' (8) 

contain consequently the hole creators bj. To conclude, in the DCM and BDCM 
one starts with the valence fermions/bosonic states of Eqs. (0, ^) and constructs 
subsequently mixed-mode nuclear states by including components having additional 
bosons formed by the particle-hole pairs of the core excitations. The resulting nuclear 
states are then classified in terms of (CMWFs) of increasing degrees of complexity 
(number of particle- hole pairs) Ref. [2]. The basic dynamic equations of the models 
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are the commutator equations that involve the nuclear many-body Hamiltonian H and 
the operators A\^{n = 1,2,3). After a lengthy but elementary algebra, we obtain the 
following results: 

(a) Commutator equation for {n = 1, 2, 3} states: 

[H,AlianJn;J)]= E {AnianJn;J)\\H\\AiiPr.Jn;J))Ali(3nJi,J) 
+ E ^ {Mc^nJn: J)||i/||4+l'(/5n+l' J)) 

xA^^^,,{/3n+iJ'JUi'^J). (9) 
(b) and Commutator equation for {n = 1, 2, 3; n' = 1'} states 

[H, A'l^-^^/^an+l'JnJn+l'] J)] 

= E {^n+l'i0in+l'JnJn+l']J)\\H\\Al^^-^,{(3n+l'JnJn+l'i'J)) 

+ E {An+v{an+vJnJn+V]J)\\H\\Al^^2'{(^n+2'JnJn+vJ'n+2')'-.J)) 

^^i+2'iPn+2'JnJ'n+l'Jn+2'i J)- (10) 

The commutator equations are used to obtain the solutions of DCM with use of 
the Equations of Motion (EOM) method [6]. In this latter method, one looks for an 
operator such that C^|0) = |m) with H\m) = Em\m) and H\0) = Eo\0) {Eo = 0). 
Here \m) denotes the excited state and |0) the correlated ground state. One then has 
the operator identity [H, C^] = EmC}^. Upon replacing C"^ by A^ and using the results 
in Eqs. |10D for the l.h.s. of the above operator identity, one obtains a set of equations 
which, after linearization, can be transformed into a system of eigenvalue equations. As 
one can see, the linearization consists in approximating the higher-order diagrams by 
an effective term. The solutions of the linearized commutator equations, can, therefore, 
be regarded as eigenvalues of a model fermionic/bosonic Hamiltonian. 

The valence particles become, therefore, the dressed solutions of the following 
collective nonlinear Hamiltonian [7]: 

N 

Hcoii = {T,oii + V,oii) = E Km AliaiJi ...Ji, J) A(«iJi ...Ji] J) (11) 

i 

with 

A\{ai{Ji . . . Ji); J) = [xl^j_jA\{aiJi ... J,; J) 

+ -^Q,+i/J,...J,+i' A>1'("*+1''^1 • • • Ji+l''-> J) + ^a^+2'J^■■■J^+2'^\+2'i^i+•i'Ji ■ ■ ■ J) (12) 

where the are the calculated mixed-mode amplitudes. 

The input to solving the commutator equations are the matrix elements of the 
n — particle configuration mixing wave functions (CMWFs). The latter can be easily 
calculated by use of the cluster-factorization method of Ref. [2, 3]. The method is guided 
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by the observation that the EOM connect the n— particle states to the {n — l)-particle 
states. In the DCM, the n ^ n — 1 reduction is achieved by factorizing out 1 boson, 
which can be {2p}, {Ip — Ih}, or {2h}. In the following, we exemplify the method 
by considering the {4p — 2h} {n = 3) parent configuration. For the sake of simpler 
notation, we will not write the detailed coupling of the creation operators but write 
the wave functions There are three active pairs associated with the $jm(o^3</i'/2'^3) 
configuration. To go to the {n = 2) (two active pairs) configuration, we can factorize out 
either {Ip— Ih}, {2p}, or {2h} to arrive at {3p— Ih}, {2p — 2h}, or {Ap}, respectively. 
This leads to the following expansion in terms of {3p — l/i}-CMWFs with coordinates 
{02}, {2p — 2/i}-CMWFs with coordinates {A2}, and {4p}-CMWFs with coordinates 
{^2}: 

\<l>T'{aMJ3)) = 

= ^ _ E 'rj(a3JlJ2J3|}«2J.«2Js) [|^S''(a2Jfc,JfcJ) ® |^jf'(«2J.))''' 

+ ^ _ E Vj(a3JlJ2J3|}A2Jj2Js)[|^r'(^2JfciJfcJ) ®|^j!(A2J.)) ' 

+ 4^ E %(a3JiJ2J3|}e2Jr.e2J.) [|^t(e24,4j) ® |^'j!:(e2Js))]'' . (13) 

V t2<^2Jk^Jk2JrJs 

The cluster transformation coefficients in Eq. (pTSD, i.e: the ^Tj{azJiJ2J-i\}ci2Jr'oi2Js)^ 
the ^Vj(a3 Ji J2<-^3|}A2^rA2 Js), and the ^ Zj{a^JiJ2Jz\}e2Jr'^iJs) are then calculated, as 
in Ref. [2, 3], by reducing the 5*^27+1(3) representations carried by the wave functions 
on the right-hand side of Eq. (^, z.e., by diagonalizing the Casimir operators [8] 



of the 5'[/2j+i(3) groups in the basis states of Eq. (p!3D . Hence, the coefficients 
^T,]{azJiJ2Jz\}<^2JrCi2Js) eigenvalues of the following matrix: 

(^Tj{a:iJiJ2h\]a2Jra2Js)f ^^/(aa Ji J2^3|}a2^ra2^s) = 

^_^yk + .J[ + J^.+jf + M, (^)l/2 Jj,^ 




(^^Tj^{a2JkiJk2\}<^lJi<^lJr)) '^Tj^{a2JkiJk2\}(^lJi<^l'Jr 
{'T\{Mk:M}(3lJ'Ajl)y 'T%{P2JkM}(3lJ'AJr) ■ (14) 



In Eq. (0), the sum over J^^, J^j, J^g, J^^, ai, oT, and jSi has not been written 
explicitly but is understood. The 6 J coefficients are defined according to Ref. [9]. The 
same procedure has been applied to calculate the and cluster-transformation 
coefficients. The cluster-transformation coefficients of the '^T are then expanded in 
terms of using Eq. (|g) in the limit -^^ T and t. In tables to | the 

transformation coefficients are exemplified. The above recursive procedure is equally 
valid for cluster transformation coefficients involving fermionic CMWFs. Within the 
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cluster transformation coefficients introduced for the fermionic/bosonic CMWFs, the 
matrix elements needed for the EOM can be readily obtained: 
The matrix elements for two particles are 

{a){^J\V\^J) = {ef\G{E)\^J) 

{h){'^'r\V\'^T) = {'^T'lV.Henl^T') (15) 

In Eq. ([l5|), G(E) is the realistic G-matrix, and Vphen is a phenomenological potential. 
The matrix elements for: 

(c){3p}, {2p - Ih}, {2h - Ih}, {3p - Ih}, {Ap}, {Ap - Ih}, {Ap - 2h} (16) 
are given symbolically in Fig. ^ 

3. HFS in hydrogen- and lithium-like ions 

With the reduction of the CMWFs and the factorization of the n-body matrix elements 
at our disposal we have performed calculations for the hyperfine splitting (HFS) of heavy 
ions [10]. In table ^ the DCM results are compared to ground-state HFS splittings 
calculated for a point nucleus and with results which take the finite spatial distribution 
of the nuclear charge into account (Breit-Schawlow correction). Additionally, QED 
radiative corrections [15] are listed and combined with the DCM results. The pure 
DCM splittings agree remarkably well with the experimental values, which are given 
in the last row of the table, while the wavelengths obtained after adding the QED 
contributions show a systematic shift to larger wavelengths. To clarify this still open 
point and in order to obtain further informations about the hyperfine interaction of the 
nucleus with the electronic cloud of high-Z ions , we have also performed calculations for 
the HFS of the lithium-like ions. The result of the calculation [16] of the HFS for 209gj80 
is given and compared with other theories in table ^. In table |^ the experimental result 
of Ref. [19] is also given. A new experiment will be performed at GSI (2003) in order 
to localize the resonance within an even smaller error. The boiling of the QED vacuum 
terms given in table 0, i.e. the electron-positrons contributions up to now not considered 
by other theories, may of course also generate a double counting effect with the QED 
contributions in lithium-like ions [16]. In the DCM the mixed-mode (b) states generate 
contributions to the HFS that may be responsible for a double counting effect with the 
QED perturbative calculation. The terms (b) are not considered by any perturbative 
calculation performed for the magnetic structure of nuclei. 

4. Elastic proton scattering on exotic nuclei 

The transformation coefficients find also their applications in the calculation of the 
matter distribution of halo nuclei [20]. Here the calculation are performed for ^Li 
and ^^Li considering a large dimensional space and treating consistently the interaction 
of the valence particles with the core excitation. For ^^Li, in order to reproduce the 
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three-body force we have performed calculations with three dressec? neutrons [21], which 
are solutions of the symbolic eigenvalue equation given in Fig. ^. The results of the 
calculations are given in Fig. ^ and Fig. ^. In Fig. ^, we plot: (left) the experimental 
matter distribution calculated in Ref. [3] with the BDCM and the experimental matter 
distribution; and (right) the ^^Li experimental scattering cross sections [22] compared 
with the results obtained in Ref. [3] for the theoretical cross section. In Fig. ^ we 
plot: (left) the calculated [3] and the experimental matter distributions; (right) the 
calculated ^^^Li scattering cross sections [21] compared with the experimental one [22]. 
For ^Li the calculated and the experimental matter radii given in Fig. |^ are in good 
agreement with the experimental matter radius given in the figure. The calculated 
matter radius of ^^Li is 3.64 fm [21] also in good agreement with the experimental value 
of 3.62(19) fm [22]. The BDCM and the DCM results obtained for the scattering cross 
sections are in good agreement with the data of Ref. [22]. For ^^Li the agreement is 
better for momentum transfer smaller than .04 (GeV/c)^, as for momentum transfer 
between 0.04-0.05 (GeV/c)^. This is probably due to high-lying core states, which are 
not included in the present DCM calculations, but are of increasing importance at higher 
momentum transfer. The consideration of these core states will increase drastically the 
dimension of the CMWFs; in the calculations presented here already about two hundred 
components were used. However, such calculations presently under investigation become 
easily feasiable with the use of cluster transformation coefficients. 



5. Proton polarization 

The scattering cross section of low-energy photons by a structureless charged system is 
given by the Thomson formula which derives the photon scattering cross section only in 
terms of linear contributions which are proportional to the frequency of the incoming 
photon. A second class of terms which are linear in the frequency of the incident photon 
and take care of the magnetic moment and spin of the target have also been included 
as corrective factor to the Thomson scattering. Additionally to these linear effects, 
however, the cross section get contributions from terms proportional to the square of 
the photon frequency. These terms, are the polarization terms of Ref. [23] and modify 
the cross section accordingly to the square bracket of the following equation: 



Ipoint ~l~ 



— — 1 + cost/ H — 1 — cos( 



E, J {hcf V 2 ' ' ' 2 ^ ' )\ ^ ' 

In Eq. (|13), the first term {^)point give the scattering of photons from a point-like 
hadron corrected with the size-spin contributions, r is the classical radius of the proton, 
Ei is the incidental photon energy, Eg^ the photon scattering energy, 9 the scattering 
angle, a the electric- and {3 the magnetic-polarizability i.e.: terms quadratic in the 
photon frequency. In Fig. |^, we plot the polarization cross section using for a and (3 the 
values calculated with the model independent dispersion sum rule in Ref. [24]: 

[a + (3) = (14.2 ± 0.03) x lO'^^m^ (18) 
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Table 1. CocfBcients for the reduction ip => {pp)j^ ® {p)js 

^3,3((3)Ji(3),Ti(l);(5,5)J2(l),T2(0)|} (5,5)J,= 1,T, = 1; (3) J, = 3,T, = 1) -0.1362 

(5, 5)Jr = 3, Tr = 1; (3) J, = 3, = 1) -0.0422 

(5, 3) J,. = l,Tr = 1; (5) J, = 5,r, = 1) 0.8613 

(5, 3)Jr = 2,Tr = 1; (5) J, = 5, = 1) -0.4876 



Table 2. Coefficients for the reduction 3plh {ph)j^ ® {pp)js 
Ti,o((3,5)Ji(2),Ti(0); (1, 1) J2(l), r2(0)|} (3,1) J, = l,r, = 0; (5, 3) J, = 1,T, = 0) 0.9326 

(3, l)Jr = 2,Tr = 0; (5, 3) J, = 2, T, = 0) 0.3606 



Table 3. Coefficients for the reduction 3plh {pp)jr ® {ph).j, 
yi,o((3, 5) Ji(2), Ti(0); (1, 1) J2(l), T2(0)|} (3, 5) J, = 1, = 1; (3, 1) J, = 2, = 0) 0.1397 

(3, b)Jr = 3, = 1; (3, 1) J, = 2, = 1) 0.9902 



These values are around those measured for the proton electric- and magnetic- 
polarization in Ref. [25]. The plotted polarization terms are ca. twelve order of 
magnitude smaller then the Compton term. Under this consideration with the present 
parameters of the PHELIX-laser in mind it is difficult to perform measurements to 
obtain the proton polarizability within a better resolution. However, additionally to the 
polarization effect produced by the energy of the incoming photon field, an intrinsic 
polarization effect [26] may contribute to the scattering cross section. This would 
result from effects analogous to the boiling of the QED terms. The three quarks are 
exciting from the vacuum in a nonlinear mechanism the {qq) condensate. The three 
dressed quarks are also eigenvalues of the equation symbolically represented in Fig. 
Calculations are in progress using for the quark the harmonic oscillator base proposed 
in Ref. [27]. The QCD boiling terms should then give rise to a radial distribution 
with an halo shape as calculated in Fig. § for ^^Li. De facto the intrinsic polarization 
effect increases the polarizability of the proton and may produce a sizable effect for 
the present parameters of PHELIX-laser. Therefore, considering this new source of 
polarization, experiments can reveal unexpected results. 

6. Conclusion 

The linearization approximation and the cluster transformation coefficients are 
important tools in solving nonlinear systems such as those characterizing n dressed 
particles: the linearization ansatz defines a collective hamiltonian that is selfconsistently 
solvable for the model eigenvalues, while the transformation coefficients provide easy 
computation for the complex n-body matrix elements which are the input to the 
collective eigenvalue equation. Within the linearization ansatz and the transformation 
coefficients, new theoretical insight have been obtained for the physics of interacting 
particles which coexist with the excitations of the model vacuum. 
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Table 4. CocfBcients for the reduction Aplh =^ {ip).j^ 'S) {ph)j, 
y3,3((l);(l,l)Ji(l),Ti(0);(5,l)J2(2),r2(l); J3(1),T3(1)|} 
(1,1,1) J,. = l,Tr = 1; (5,1) J, = l,Ts = 1) 0.9513 
(1,1, 5) J,. = 3,Tr = 2; (1,1) J, = l,Ts = 1) 0.3081 



Table 5. CoefRcients for the reduction Aplh ^ {2plh)j^ ® {pp)j, 
iy3,3((l)(3,5)Ji(l),Ti(0);(5,l)J2(2),T2(l); J3(1),T3(1)|} 
(1,1,1) J, = l,r, = 1;(1,5)J, = 1,T, = 1) 1.0 



associated 



associated 



associated 



configuration matrix elements / potentials configuration matrix elements / potentials configuration matrix elements / potentials 



2p |J 
7p-1/l 1^ 



realistic G matrix 



phenomenological particle- 
iiole matrix element 



realistic G matrix 



3p ^ 

=z kV\\ 



ijk i j k 

realistic G matrix and 



2p-1h phenomenological 



particle-hole potential 



ijk 



4 



I J k 



i jk 



4p 



realistic G matrix 



-'ijkl 

ijkl i j k I 

I i realistic G matrix and 
3p-\n phenomenological 
r I particle-hole potential 

ijkl |_ ijkl ijkl 

I i I realistic G matrix and 
4p-1n phenomenological 
f^'^ particle-hole potential 



ijkimn 



ijkim 



Qjklr^^ Q^m 

I j k I m i i k Im _ 
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Figure 2. Simbohc representation of dressed three particles. 
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Figure 4. Left: Experimental and calculated distributions of ^^Li; Right: 
experimental and theoretical proton scattering cross sections. 
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Figure 5. Scattering cross sections of high energy photons as functions of scattering 
angle and proton polarization parameter. 
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Table 6. Ground state hyperfine structure splittings for a point-nucleus (-Bpn), 
including Breit-Schawlow (Ebs), DCM, and QED corrections. The QED contributions 
include vacuum polarization and self energy. All values are in eV. 





185j^g74+ 


187j^g74+ 


203rp|80+ 


205rp|80+ 


207p]-,81+ 


209gj82+ 




3.0103 


3.0411 


3.0184 


2.9890 


1.3998 


5.8395 


Ebs 


2.7976 


2.8263 


3.3073 


3.3374 


1.2528 


5.1922 


TpDCM 


2.7192 


2.7449 


3.2130 


3.2770 


1.2166 


5.0832 


A-Eqed 


-0.0142 


-0.0143 


-0.0176 


-0.0177 


-0.0067 


-0.0280 


Etot 


2.7050 


2.7306 


3.1954 


3.2213 


1.2099 


5.0552 


E 


2.7187 (18) 


2.7449 (18) 


3.21351(25) 


3.24409(29) 


1.2159 (2) 


5.0841 (4) 


Ref. 


[11] 


[11] 


[12] 


[12] 


[13] 


[14] 



Table 7. Comparison of different hfs calculations in meV for Li- like bismuth with 
experimental HFS splitting 



Contribution (meV) 


[17] 


[18] 


[16] 


one-electron 


958.50 (5) 


958.49 


958.51 


charge distr. 


-113.8(3) 


-151.44 


-113.61 


mag. distr. 


-13.9(2) 


-18.68 


-14.1 


total QED 


-4.44 


-4.06 


-4.81 


e-e interaction 


-29.45 (4) 


8.43 


-3.4 


boiling of QED vacuum 




0.06 


-7.69 


total theory 


796.9(2) 


792.8 


783.9(3.0) 


measurement [29] 




820 (26) meV 
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